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We study the cosmological evolution of an induced gravity model with a self-interacting scalar 
field a and in the presence of matter and radiation. Such model leads to Einstein Gravity plus 
a cosmological constant as a stable attractor among homogeneous cosmologies and is therefore a 
viable dark-energy (DE) model for a wide range of scalar field initial conditions and values for its 
positive 7 coupling to the Ricci curvature ^a 2 R. 

PACS numbers: 



I. INTRODUCTION 



Several years ago a model for a varying gravitational 
coupling was introduced by Brans and Dicke The 
model consisted of a massless scalar field whose inverse 
was associated with the gravitational coupling. Such a 
field evolved dynamically in the presence of matter and 
led to cosmological predictions differing from Einstein 
Gravity (EG) in that one generally obtained a power- 
law dependence on time for the gravitational coupling. 
In order to reduce such a strong time dependence in a 
cosmological setting, while retaining the Brans-Dicke re- 
sults in the weak field limit, several years ago a simple 
model for induced gravity @, [|[ involving a scalar field 
a and a quartic Act 4 /4 potential was introduced. This 
model was globally scale invariant (that is did not in- 
clude any dimensional parameter) and the spontaneous 
breaking of scale invariance in such a context led to both 
the gravitational constant and inflation, through a non- 
zero cosmological constant The cosmological con- 
sequences of introducing matter as a perturbation were 
studied leading to a time dependence for the scalar field 
and consistent results. Since the present observational 
status is compatible with an accelerating universe domi- 
nated by dark energy (DE) we feel that the model should 
be studied in more detail. 

In an EG framework quadratic or quartic potentials for 
canonical scalar fields are consistent with DE only with 
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an extremely fine tuning in the initial conditions leading 
to slow-roll evolution until the present time. Indeed, a 
massive or a self-interacting scalar field in EG, respec- 
tively behave as dust or radiation during the coherent 
oscillation regime. In contrast with this, induced gravity 
with a self-interacting Act /4 potential has as attractor 
EG plus a cosmological constant on breaking scale in- 
variance. The simple model we consider illustrates how 
non-trivial and non perturbative dynamics can be ob- 
tained in the context of induced gravity DE, and more 
generally within scalar-tensor DE. 



II. THE ORIGINAL MODEL 



Let us consider a system described by the Lagrangian 



C=-^i,{-g^d^d v a + 1 a 2 R~-a i \+ £ C 3 

\ ' j=R,M 

(1) 

where 7, A are dimcnsionlcss, positive definite parameters 
and the C% 's are the contributions of cosmological fluids 
behaving as dust and radiation. If we consider the homo- 
geneous mode for the scalar field evolving on a spatially 
flat Robertson- Walker background 



ds 2 = g^dx^dx" = -dt 2 + a 2 (t)dx 2 (2) 
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from the above Lagrangian one obtains the following set 
of independent equations 

' H 2 = V ^3 + 7T 4 -2H2- + ^-(T 2 

<i-— ' 3-7cH 6-7 <7 2 cr 12-7 

, p i = -ZH{ Pj +P i ) 

where P 3 - = WjPj, wr = 1/3, ium = 0, the dot de- 
notes differentiation with respect to the cosmic time t 
and H = a/ a. Note that the system ([3]) of equations are 
respectively the Friedmann, the Klein-Gordon and the 
continuity equations. In the above context, the critical 
densities for radiation, matter and a can be taken to be: 

O PR O PM 

(4) 

- _ (j 2 -12gg-o-+^) 
— (6~/<r 2 H 2 ) ' 

where f2j\/ + + CT = 1 . 

On using the equations of motion one finds that the cur- 
vature scalar is given by 

W 2 y 1 + 67 J 

Even in the absence of the fluids pj , the quantity {5J has 
a undefined sign, implying that the a potential has one 
or two minima depending the phase space trajectory of 
a. 

We also observe that the above model ([TJ) corresponds to 
an interacting Brans-Dicke field 4>. Indeed on redefining 
cf> = a 2 one has 
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FIG. 1: Phase space plot of equation Q); on the horizontal 
axis the field <p is represented while on the vertical we plot 
its derivative 4j& with 7 = 2- 10 _ , A = 10 _1 , co = 1 and 
r ( ) = (dotted trajectories) or t~o = 5 • 10" 1 (dashed trajec- 
tories). The solid lines are the boundaries of the accessible 
region of the phase space while the dotted-dashed line repre- 
sents all the possible solutions with a = 0. The shaded area is 
the portion of phase space where the scalar field has a double 
well potential. 



The possibility of super-acceleration is due to the ^Haa 
contribution in the brackets. Let us note that 7 should be 
large enough in order to have a super-accelerated phase 
which can be distinguished from a de Sitter phase. This 
cannot be accomplished in EG with a scalar field having 
a standard kinetic term. 



B. Scalar Field plus Radiation 



where oj — (47) -1 . The above model {1} also admits a 
transformation in the Einstein frame in which the po- 
tential for the redefined field is just a cosmological con- 
stant (= 167r 2 G 2 A/7 2 ) and dark matter and the redefined 
(free) scalar field are coupled (by (7 2 ct 4 ) -1 ). 



A. Super-acceleration 

Let us note that scalar-tensor theories admit naturally 
super-acceleration, i.e. H > (see Q for a generalization 
of the model in 0,H[). From the Friedmann equation and 
expression ([5]) one finds 



H = — 



1 

ja 2 



E 

j=R,M 



(l + 87 + Wj ) 
2(1 + 67) 



-Pj 



cr 
~2 



A'yHcrcr 



•(7) 



When pressureless dust can be neglected with respect 
to radiation (pn 3> pm) & phase space analysis can be 
performed for the system. In this case, the Klein-Gordon 
equation in ([3]) can be integrated and leads to 



H 



(8) 



where a prime denotes a derivative with respect to N = 
In a and Co is an integration constant. The Friedmann- 
like equation can then be cast in the form of an au- 
tonomous equation describing the possible solutions for 
the rescaled field ip = a a: 



<P 



1 (^-^) 2 + 2^(^-^) 



6-7 



(9) 
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where we set Pr$ = ArL=i (similarly we shall set 
Pm ,o = /°A/| a =i)- The right hand side of eq. d9j) is a 
semi-positive definite quantity and thus its left hand side 
must satisfy the inequality 

ip 2f3 = ip *-L( (p >- (p f+2<p(<p'-<p)>0 (10) 

which determines the allowed region of phase space. 
Once 7 and A are fixed, one must choose the initial condi- 
tions encoded in cq and pr.q in order to determine com- 
pletely a phase space trajectory. Note that if 7 = 
or Co = there is only one possible trajectory given by 
(p' = (p corresponding to the GR limit a = 0. 



III. NUMERICAL ANALYSIS 
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FIG. 2: Plots of , Qm , Ode for the sets of parameters and 
initial conditions 1 and 2 of Table I. 



If a realistic cosmological evolution is considered with- 
out any approximations, the system ([3|) can only be 
solved numerically. On defining a rescaled Hubble pa- 
rameter TL = \~ 1 / i H and introducing the dynamical 
variable x = A 1 / 2 ^ 2 the relevant equations can be rewrit- 
ten as 



H 2 x" 



H 



1 - 



2H 2 + \(H 2 )' 



67+I 



_ Pm+Pr 1 x 
37X 127 



(11) 



with pr — pRfl/a A and pu — Pm.o/o. 3 - In this form 
the A dependence simplifies and one can find numerical 
solutions with a high precision even when A is extremely 
small which is the case (note that this is crucial since 
in the A — > limit one would otherwise just recover the 
Brans-Dickc theory). 

The initial kinetic energy-density of the a field is rapidly 
dissipated during early stages when radiation dominates 
over dust, as can be deduced from phase space analysis 
and is confirmed by Fig. {1}. Consideration of high val- 
ues for the time derivative of the scalar field (when it 
dominates the energy-density of the universe, see Fig. [2|) 
lead to a "stiff" universe. If it existed, such an epoch 
should have occurred well before nucleosynthesis in or- 
der not to spoil the observed abundances of the light 
elements. Indeed this could be the case for very high 
densities arising from fermions (quarks) interacting with 
vector mesons (gauge fields) 

The value for a remains constant until the onset of the 
matter dominated era and then its variation is well ap- 
proximated by 



rfo 2 
dN 



I + 67 



(12) 



where we use cr = c| z=0 , and then it approaches the 
constant value finally leading to the dc Sitter stage. The 
above relation (|12p is obtained by assuming a matter 
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FIG. 3: Evolution of (3 defined in U0\) plotted for the sets of 
initial conditions and parameters 1 and 2 of Table I. Note that 
the solid and the dotted curves are indistiguishable between 
N 8 and N ~ 



(dust) dominated universe and is as expected within the 
context of induced gravity theories [l], Q . 

From here on, we shall consider different choices of pa- 
rameters and initial conditions for our analysis (see table 
(|T])). In the first two cases in table (JXJ) 7 has been taken to 
be its maximum value allowed by solar system constraints 
[f| 0, [H, while initial conditions for a describe two ex- 
treme situations: fii n i,i — and (3mi.2 — 1- The early 
stages of the evolution, when the matter contribution is 
negligible, are well approximated by the plots in figure 
(TTJ. As expected, the field a rapidly evolves toward the 
shaded region where j3 ~ 1 (see Fig. (J3j) ) . Since 7 << 1, 
when j3 ~ the Hubble parameter scales as a~ 3 . 
Let us end this section by noting that one should not 
evolve the scalar field equations backwards in time - i.e. 
with rcdshift - since the correct physical trajectory on 
which the scalar field evolves can not be identified. 
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TABLE I: Parameters, initial conditions and results for different cases. All these cases have Qm ~ 0.27 , Ade — 0.73, Ho = 
73 kms _1 Mpc~ 1 , 7a"o — (SivG)^ 1 ,jm = 5 • 10~ 7 . Note that the solar system bound on 7 is relaxed for # 3 to 8; nonetheless 
# 3 to 7 of the model are still satisfactory on comparing with cosmological data coming from SNIa and first Doppler peak of 
CMBR anisotropies. The smaller the value of 7 the closer is our model to ACDM. 



A. Solar system Newtonian constraints 

The parameters of a scalar-tensor theory can be con- 
strained though solar system observations which test 
classical general relativity [§] . Since our model coincides 
with the Brans-Dicke case in the weak field limit 0, Q , 
one has [1, [T3| Ppn — 1 = and 



1PN - 1 = - 



47 



I + 87 



(13) 



Therefore, since our 7 is positive and "/pn — 1 = (2.1 ± 

2.3) -lo- 5 we shall take an upper bound for 7 given 

by 7m = 5 ■ 10~ 7 which we shall use. Let us note however 
that for other values of -fpN reported in the body of Ref. 
[(| the resulting 7m would be 5 to 10 times larger. 



In scalar-tensor gravity the effective Newtonian constant, 
namely the gravitational coupling measured in Cavendish 
like experiments, is 



G, 



1 



ff,0 



87 + I 



8tt 7 Cg 67 + 1 



(14) 



Possible variations of the effective Newtonian constant 
G e e,o, are allowed in the interval —7 • 10 _14 y _1 < 
G cff , /G cff ,o < 3 • 10~ 14 j/- 1 0. The value G N of New- 
ton's constant is known with a precision of ~ 10 -3 [Til ] ; 
we shall thus restrict our analysis to a set of initial con- 
ditions for a leading to values of G e ff,o within the bounds 
given by the above experimental uncertainty on the mea- 
surement of Gn- 
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IV. COMPARISON WITH EINSTEIN GRAVITY 



In order to compare a scalar-tensor theory with 
dynamical scalar field models in EG (also dubbed 
quintessence [12]), we first identify DE energy-density 
and pressure. A consistent definition of such quantities 
can be given [l(| if one identifies the effective energy- 
density and pressure in the set of equations (|3|) in an EG 
framework with 8itGn = 7<7g: 



' Pde = ZlolH 2 - Pj 

j=R,M 

Pde + Pde = - J2 (pj + Pj) 

j=R,M 



(15) 



Such an identification mantains the continuity equation 
for DE and leads explicitly to 
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■2-YPj+Pj 
I+67 



PR 



(16) 

We can as well identify the relative densities for the fluids 
involved in the evolution 



PR 



3 7 <7 2 iP 



, n 



M 



PM 



37^# 2 ' 



Ode = 



Pde 
3 7( t 2 # 2 



(17) 



where Qm + ^R + ^de = 1- We stress that these crit- 
ical densities based on this fictitious EG model differ 
from the previous ones ^ wherein the gravitational cou- 
pling varies. For 7^1 differences between the two 
definitions are negligible. The matter and DE critical 
densities are plotted in Fig. The state parame- 

ter wde defined in Eq. p^|) is displayed in Fig. §5§ 
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for 7 = 5 x 10" 7 , 1CT 4 , 10" 2 (cases # 1, 5, 8 in Table 
(jl])). Note that interpreting this self-interacting induced 
gravity as an EG model, the parameter of state of DE 
is 1/3 during the radiation era, is positive during the 
matter dominated stage (with a non negligible bump for 
7 = 1CP 2 ) before setting to —1 + 0(7) at present, as can 
be seen in Table I, where we have reported few examples 
of the simulations we have performed. Let us also note 
how the EG model is not completely dominated by mat- 
ter, in particular for 7 ~ 10" 2 . 

Since the state parameter at present differs only slightly 
from —1 it is therefore appropriate to compare this self- 
interacting induced gravity model with a ACDM model 
with f^DE .0 a s the fractional energy-density stored in the 
cosmological constant. The relative deviation of such an 
Einstein ACDM model from the original model for the 
comoving distance 

d L {z) = (i + z) J* (18) 

evaluated at the last scattering surface (z = 1089) is 
shown in the last column of Tab. (jl|. The comoving dis- 
tance to the last scattering surface enters in the acoustic 
scale, which is related to the characteristic angular scale 
of the peaks of the CMB angular power spectrum. The 
current uncertainty on the acoustic scale is of the order 
of one percent; on the basis of the numbers reported in 
Table QJ, this implies that models with 7 < 0(1O~ 2 ) 
may be constrained by the position of the peaks of the 
CMB angular power spectrum. 

We now compare the model with data from Supernovac 
la, which constrain the evolution of the scale factor for 
< z < 2, as done in other works studying scalar-tensor 
DE [l3|. The supernova data we use are the "gold" set 
of 182 SNe la by Riess ct al. [l~|. In scalar-tensor theo- 
ries, since the Chandrasckhar mass (~ x 3 / 2 ) varies with 
redshift, the modulus-redshift relation is: 

H(z) = 5 log 10 — — + — log 10 — - . 19 
10 pc 4 x(z) 

In table Q the x 2 values are listed showing that 7 ~ 10~ 2 
is strongly disfavoured by SNe la. An accurate statistical 
analysis is clearly beyond the scope of this work, but we 
expect that SNe la may constrain 7 ~ O(10~ 3 ). 

V. CONCLUSIONS 

We have shown how a simple model of self-interacting 
induced gravity @, [|[ is a viable DE model, for tiny 
values of the self-coupling (A ~ O(10" 128 )). The model 
has a stable attractor towards EG plus A and can be very 
similar to the ACDM model for the homogeneous mode, 
on taking into account the Solar system constraints 
quoted in [1, 0, 0- At the cosmological level, in the 
presence of radiation and dust, it is interesting that 
for such a simple potential (quartic for induced gravity 
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FIG. 4: Evolution of Om and Qde as defined in (|16[) for dif- 
ferent choices of 7. 
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FIG. 5: Evolution of wde as defined in (|16|l for different 
choices of 7. 



or interacting for the equivalent Brans-Dicke model) 
the model has an attractor towards EG plus A, very 
differently from the case of a massless scalar - i.e. A = 
-, for which there is no mechanism of attraction towards 
EG. The attractor mechanism towards GR and the 
onset to acceleration are both inevitably triggered by 
the same mechanism, i.e. scale symmetry breaking in 
this induced gravity model. In contrast with EG, the 
choice of a runaway potential Il5j for quintessence is 
not mandatory (see however [16[ for a study of these 
potential in the context of scalar-tensor theories). If 
the final attractor is an accelerated universe, the only 
constraints on parameters and initial conditions come 
from observations. We have shown that late time 
cosmology - after recombination, for instance - is mostly 
insensitive with respect to the initial time derivative of 
a. For this reason, the full set of parameters and initial 
conditions of the model are fully specified on taking the 
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observed values for G, Hq, 

We have discussed in detail such model in the context 
of Einstein gravity, i.e. keeping the Newton constant 
(approximately) fixed at its actual value. The model 
predicts an equation of state of the equivalent Einstein 
model with wde slightly less than —1 at present and 
homogeneous cosmology by itself seems able to constrain 
7 < 0(1O~ 3 ), although a full statistical analysis is 



clearly beyond the scope of this work. 
It is interesting to also study also what other potentials, 
besides the simple potential Aer 4 /4 employed in this pa- 
per, will also be compatible with the observed evolution 
of the universe. 
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